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The  progress  made  during  the  third  year  of  an  investigation  to 
measure  the  response  of  a burning  solid  propellant  to  oscillatory 
flow  conditions  is  presented.  In  this  study  a modification  of  the 
impedance  tube  technique  is  used  to  measure  the  response  over  a 
wide  range  of  frequencies.  Improvements  in  the  data  reduction  pro- 
gram are  discussed.  These  include  a more  accurate  method  of  com- 
puting the  temperature  distribution  in  the  burner  tube  and  a tech- 
nique for  determining  the  response  from  pressure  amplitude  measure- 
ments only.  A high-pressure  facility  and  minicomputer-based  data 
acquisition  system  are  also  discussed  in  the  report.  Data  taken  at 
300  psig  indicate  that  the  combustion  process  of  the  solid  propel- 
lant periodically  drives  and  damps  acoustic  oscillations  under  most 
of  the  test  conditions  encountered. 
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I.  INTRODUCTION 


This  report  summarizes  the  work  done  during  the  third  year  of 
research  supported  by  Air  Force  Grant  No.  AFOSR-73-2571 , initiated 
July  1,  1973.  This  grant  is  concerned  with  the  measurement  of  the 
acoustic  responses  of  solid  propellants  by  the  impedance  tube 
me thod. 

The  acoustic  response  is  an  important  parameter  in  the  study  of 
combustion  instability  in  solid  propellant  rockets.  Combustion  in- 
stability is  characterized  by  high  amplitude  pressure  oscillations 
which  are  driven  by  the  combustion  process.  These  high  amplitude 
pressure  oscillations  can  cause  large  variations  in  the  thrust  with 
time;  can  interfere  with  delicate  control  and  guidance  systems;  and, 
in  extreme  cases,  can  lead  to  structural  failures.  The  acoustic  re- 
sponse is  a measure  of  the  capability  of  a specific  burning  solid 
propellant  to  initiate  and  sustain  combustion  instability  in  a given 
solid  propellant  rocket  motor. 

To  analyze  the  stability  of  a particular  motor,  it  is  necessary 
to  determine  the  response  (or  reaction)  of  the  combustion  process  to 
the  various  disturbances  (e.g.,  pressure  fluctuations  associated  with 
turbulence  or  combustion  noise)  inside  the  combustion  chamber.  The 
amplification  or  decay  of  these  fluctuations  will  depend  upon  the 
phase  relationship  between  the  pressure  oscillations  and  the  induced 
combustion  oscillations. ^ If  the  response  time  of  the  combustion  pro- 
cess to  the  combustion  pressure  oscillation  is  such  that  an  increase 
in  burn  rate  occurs  during  a pressure  condensation,  then  a condition 
supporting  an  increase  in  the  amplitude  of  the  pressure  oscillation 


1 


1 


may  result,  possibly  leading  to  an  unstable  engine  operation.  The 


l l 

m 

| i 


question  whether  instability  will  occur  also  depends  upon  the 
"losses"  in  disturbance  energy  that  are  present  in  the  combustor. 

On  the  other  hand,  if  the  burn  rate  increase  occurs  during  a pressure 
rarefaction,  a condition  leading  to  destructive  interference  will 
result  and  the  disturbance  will  be  damped. 

The  combustion  response  can  be  related  to  the  admittance  of  the 
burning  surface  which,  in  turn,  can  be  used  to  evaluate  the  effect 
of  the  burning  propellant  upon  the  stability  of  the  solid  rocket 
under  consideration.  Tiffined  as  the  ratio  of  the  normal  velocity 
perturbation  to  the  pressure  perturbation  at  the  burning  propellant 
surface;  a positive  admittance  indicates  that  the  normal  velocity 
and  pressure  oscillations  are  in  phase  and  a situation  leading  to 
disturbance  amplification  may  result.  Negative  admittance  values 
indicate  that  the  normal  velocity  and  pressure  perturbations  are 
out  of  phase  and  damping  of  the  disturbance  at  the  burning  propel- 
lant occurs.  Once  determined,  the  admittance  can  be  used  as  a bound- 
ary condition  in  combustion  instability  analysis  to  predict  the  sta- 
bility of  a given  engine  configuration. 

The  most  widely  used  method  for  measurement  of  the  admittance 

is  the  T-burner.  However,  several  difficulties  with  this  method 

have  been  encountered  and  are  discussed  in  Ref.  2.  The  present  study 

is  concerned  with  the  development  of  an  impedance  tube  method  for 

3 

measuring  the  admittances  of  burning  solid  propellants.  This  tech- 
nique, described  in  Ref.  2,  is  known  to  yield  accurate  results  in 

4 

studies  associated  with  sound  absorbing  materials  and  has  been 
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successfully  used  in  the  study  of  the  damping  provided  by  typical 
solid  rocket  nozzle  configurations.^ 

In  this  report  progress  during  the  past  year  is  discussed  and 
improvements  in  the  supporting  theory  of  the  impedance  tube  technique 
are  presented.  These  include  a more  accurate  method  for  determining  the 
steady-state  temperature  profile  in  the  tube  and  a method  for  com- 
puting the  admittance  from  pressure  amplitude  measurements  only.  In 
addition,  the  development  of  a high  pressure  facility  capable  of  meas- 
uring admittances  at  pressures  up  to  500  psig  is  described.  The 
computerized  data  acquisition  facility  which  became  operational  this 
past  year  is  also  discussed  and  preliminary  admittance  results  at 
300  psig  are  presented. 
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ANALYTICAL  INVESTIGATIONS 


The  equations  describing  the  behavior  of  the  standing  wave  pattern 


inside  a driven  burner  tube  are  discussed  in  detail  in  Ref.  3.  Also 


presented  are  the  techniques  used  to  obtain  the  admittance  of  a burn 


ing  solid  propellant  from  experimentally  measured  pressure  amplitude 


and  phase  data  taken  at  discrete  points  along  the  standing  wave.  To 


further  improve  the  accuracy  of  the  computed  admittance  values,  addi 


tional  refinements  have  been  incorporated  in  the  data  reduction  meth 


ods  during  the  reporting  period,  and  these  are  described  in  this  sec 


tion 


A.  Determination  of  the  Steady  State  Temperature  Profile 


In  Ref.  3 a method  for  determining  the  steady  state  temperature 


distribution  in  the  driven  burner  tube  is  presented.  From  this  temp 


erature  distribution  the  variation  in  wavelength  along  the  axis  of 


the  tube  can  be  determined  and  accounted  for  in  the  admittance  com 


putations.  The  equation  for  the  temperature  distribution  is  obtained 


from  the  steady  state  energy  equation  in  the  tube  * and  is  given 


by  the  expression 


temperature  at  the  wall  of  the  tube 


chamber  temperature 


Since  T and  T are  known  from  experimental  measurement,  the  con 


stant  C is  the  only  parameter  to  be  resolved  in  order  to  determine 


the  temperature  distribution  in  the  tube.  As  shown  in  Ref.  3,  this 


WWW  l".T» 


constant  can  be  found  by  measuring  the  distance  between  two  successive 
pressure  minima  along  the  standing  wave  pattern.  The  accuracy  of  the 
value  of  C depends  upon  how  accurately  the  positions  of  the  pressure 
minima  can  be  determined.  In  the  past  these  positions  were  determined 
from  pressure  amplitude  data  measured  by  transducers  placed  at  various 
locations  along  the  walls  of  the  burner  tube.  Since  the  locations  of 
the  transducers  in  general  did  not  coincide  with  the  locations  of  any 
of  the  standing  wave  minima,  these  minimum  positions  were  obtained  by 
extrapolation  of  the  measured  data  - a process  which  inevitably  results 
in  errors  in  locating  these  standing  wave  minima.  Hence,  a search  for 
a method  for  the  determination  of  the  steady  state  temperature  dis- 
tribution in  the  tube  was  undertaken  during  this  reporting  period. 

An  alternate  technique  which  decreases  the  inaccuracies  involved 
in  locating  the  pressure  minima  has  recently  been  developed.  This 
method  employs  the  nonlinear  regression  technique  to  determine  the 
constant  C from  the  measured  amplitude  and  phase  data.  The  nonlinear 
regression  technique  involves  determining  the  value  of  C which  pro- 
vides the  best  fit  between  the  theoretically  predicted  wave  pattern 
in  the  burner  tube  and  the  experimentally  measured  pressure  ampli- 
tude and  phase  data.  The  best  fit  is  found  by  determining  the  value 
of  C which  minimizes  the  root-mean-square  deviation  between  the 
theoretically  predicted  burner  tube  wave  structure  and  the  corres- 
ponding experimental  data.  This  technique  gives  a temperature  pro- 
file which  is  consistent  with  the  measured  wave  structure. 

To  find  the  minimum  root-mean-square  deviation,  the  following 
function  F,  which  is  a measure  of  the  error,  is  minimized; 

n ? 

F = E - T.)  (3) 

i=l 
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( 
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In  Eq.  (3),  represents  an  experimentally  measured  quantity  (i.e., 
amplitude  or  phase)  at  location  i and  is  the  corresponding  theo- 
retically predicted  quantity  which  depends  among  other  things,  upon 
the  constant  C.  If  a minimum  of  F exists  with  respect  to  C then  the 
gradient  of  F must  vanish  at  that  minimum  and  the  following  relation- 
ship holds : 


1=  ^[(vvir]-» 


di . 


(4) 


A Newton-Raphson  iterative  scheme  is  used  to  obtain  a solution  of 
Eq.  (4)  utilizing  a linearized  version  of  this  equation  which  in- 
volves the  expansion  of  in  a first  order  Taylor  series  with 

3 

respect  to  the  parameter  C.  The  resulting  linear  algebraic  equation 
can  be  expressed  in  the  following  form: 

( cm+1-  cm 


&<Ei-0-sfL  -[(  Cn+!  c”  ) £,  (3r)|2  ] 

1-1  l=1  V /l  (i,m) 


(5) 


where  m represents  the  mth  iteration. 
Equation  5 can  be  rewritten  as  follows: 

Am  [Cm+1-  Cm  ] = Bm 

Cm+1  = Cm  + Bm/Am 


(6) 


where 


and 


&(«) 


(i,m) 


(7) 
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Equation  (7)  is  a linear  algebraic  equation  for  the  unknown  C 
and  it  can  be  readily  solved  once  Am  and  Bm  are  computed.  The  compu- 
tation of  these  elements  requires  the  determination  of  the  derivative 
dT 

...  Since  these  derivatives  cannot  be  determined  analytically, 

C (i>m) 

numerical  values  for  these  derivatives  are  obtained  using  finite 
difference  methods. 

B . Computational  Scheme  to  Minimize  the  Error: 

As  discussed  in  Reference  3,  the  non-linear  regression  is  used 

to  find  the  values  of  the  velocity  wave  u7  , pressure  wave  p7  and 

o’  ro 

density  p7  at  the  propellant  surface  which  give  the  best  fit  be- 
tween the  theoretically  predicted  acoustic  pressure  distribution 
and  the  corresponding  experimentally  measured  acoustic  pressures. 
These  values  are  then  used  to  determine  the  propellant  admittance 
and  response  factor.  In  the  last  section  it  has  been  shown  that  a 
similar  non-linear  regression  technique  can  be  used  to  obtain  the 
constant  C which,  in  turn,  is  used  to  compute  the  steady  state  tem- 
perature distribution  that  is  consistent  with  the  measured  acoustic 
pressure  data.  Since  these  two  non-linear  regression  operations  are 
coupled,  they  must  be  solved  together  to  obtain  the  propellant  ad- 
mittance and  response  factor  data.  The  manner  in  which  this  is  being 
done  is  briefly  explained  as  follows: 

(1)  Initially,  an  arbitrary  value  of  C is  chosen,  and,  using  the 

3 

"Transmission  Matrix"  scheme  , a set  of  "initial"  values  for  u7  , p' 

o ’ ro 

and  p'  at  the  propellant  surface  is  computed. 


(2)  The  non-linear  regression  scheme  is  used  to  obtain  the  values 

of  u7  , p'  and  p 7 which  minimize  the  root-mean-square  deviation  be- 
o ’ ro  "o  n 

tween  the  theoretical  equations  and  the  experimental  values  for  the 
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chosen  value  of  C 


I i 


(3)  Once  these  are  determined,  the  value  of  C is  recomputed  using 

the  determined  "optimum"  values  of  u'  , p'  and  n'  obtained  earlier. 
r o ’ *o  Ko 

Steps  (2)  and  (3)  are  repeated  until  the  variations  in  the  values  of 

Uo  * Po  * po  ’ anc*  C a specified  tolerance. 

To  check  the  accuracy  of  this  data  reduction  scheme,  computations 

were  performed  for  a number  of  hypothetical  cases  for  which  the  exact 

values  of  C,  u'  , p'  and  0'  were  assumed  to  be  known.  Using  these 
o o Ko 

values  as  initial  conditions,  the  burner  tube  equations  were  solved 
to  determine  the  "exact"  wave  structure  in  the  hypothetical  burner 


tube;  the  latter  is  plotted  (for  a specific  case)  as  a solid  line 
in  Fig.  1.  Next,  the  points  denoted  by  "x"  on  Fig.  1 were  assumed  to 
be  measured  pressure  data  that  deviated  from  the  exact  data  due  to 
experimental  errors.  This  data  was  then  input  into  the  data  reduction 
programs  in  an  effort  to  determine  the  error  in  the  admittance  com- 
putation resulting  from  experimental  inaccuracies.  An  arbitrary  in- 
itial value  for  C was  chosen  and  the  computational  scheme  discussed 
earlier  was  used  to  compute  the  wave  structure  and  the  desired  value 
of  C.  The  computed  wave  structure  is  represented  by  the  dashed  line 
in  Fig.  1 and  the  table  below  provides  a comparison  between  the  "exact" 
and  computed  values  of  Pq  » Uq  » p'0  and  C 


Computed 

Values 

Guessed 


TABLE  I.  Ihe  Effect  of  Assumed  Experimental  Errors  on  C 


Heat  Transfer  Pressure  Pq 


Velocity  Uq 


Density 


Parameter 

C 

Amplitude 
in  db 

- 

Phase 
in  Deg. 

Amplitude 

Ft/Sec 

Phase 
in  Deg. 

Amplitude 

Slug/ft-* 

Phase 
in  Deg 

0.943 

145.4 

230.2 

1.48 

-21.4 

.612xl0'6 

229.0 

0.972 

145.8 

229.9 

1.30 

-18.2 

.647x1 O'6 

228.7 

2.0 
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Figure  1.  Comparison  of  an  Exact  Standing  Wave  Profile  and  a Computed 
Profile  with  Assumed  Experimental  Errors. 


Examination  of  the  data  shown  in  Table  I indicates  good  agreement 
between  the  assumed  and  computed  values  of  u1  , z>'  , o'  and  C even 
when  small  experimental  errors  are  present.  This  good  agreement  pro- 
vides further  confirmation  that  the  data  reduction  program  is  working 
satisfactorily. 

C.  Error  Minimization  by  Application  of  a Regression  Technique  Usin 


Pressure  Amplitudes  Onl' 


Before  incorporation  of  the  computerized  data  acquisition  system 
into  the  test  procedure,  it  was  found  that  the  magnitude  of  the  errors 
in  the  measured  phases  at  low  amplitudes  near  the  pressure  minima  were 
unacceptable.  This  necessitated  the  development  of  a scheme  to  obtain 
values  of  u^  , p^  and  at  the  propellant  surface  using  pressure 
amplitude  data,  which  could  be  measured  within  acceptable  tolerance 
near  the  pressure  minima. 

In  this  scheme,  the  starting  values  of  u^  , p^  and  are  computed 
3 

as  before  using  the  transmission  matrix  technique.  However,  in  this 

case,  the  non-linear  regression  is  performed  using  only  pressure  amp- 

3 

litude  data.  In  the  earlier  scheme  the  derivatives  of  both  the  real 
part  and  imaginary  parts  of  the  complex  pressure  are  obtained  with 
respect  to  the  real  and  imaginary  values  of  u^  , p^  and  p\  In  the 
modified  scheme  these  derivatives  are  used  to  obtain  the  derivative 
of  pressure  amplitude  with  respect  to  the  real  and  imaginary  values 
of  u^  , p^  and  Representing  the  real  and  imaginary  parts  of 

u'  , p^  and  by  2^,  i = 1,2,  ---,6  the  derivatives  of  pressure 
amplitude  p^  are  given  by 


aV(PR2  + v'l  > 1 / K 
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where  and  represent  the  real  and  imaginary  parts  of  the  oscil- 
latory pressure.  Once  these  derivatives  are  computed  the  non-linear 
regression  can  be  used  to  minimize  the  error  between  theoretically 
computed  and  the  experimentally  measured  pressure  amplitudes.  The 
development  of  this  modified  scheme  has  been  completed  and  it  is 
currently  available  for  data  reduction. 

III.  EXPERIMENTAL  APPARATUS 
A.  High  Pressure  Experimental  Facility 

Previous  experiments  using  the  driven  burner  tube  to  measure  the 
admittance  of  burning  solid  propellants  have  been  limited  to  chamber 
pressures  less  than  50  psig  because  a suitable  acoustic  driver  cap- 
able of  higher  operating  pressures  was  not  available.  In  addition, 
pressure  transducers  capable  of  withstanding  pressures  of  up  to 
500  psig  and  having  the  desired  resolution  (i.e.  0.01  psi  to  3 psi) 
were  not  available. 

Recently,  however,  a high  pressure  facility  designed  for  operating 
pressures  up  to  500  psig  has  been  developed.  This  facility  is  capable 
of  determining  the  admittances  of  burning  solid  propellant  samples 
under  test  conditions  simulating  actual  rocket  motor  operating  con- 
ditions. A photograph  of  the  high  pressure  driven  burner  facility  is 
shown  in  Fig.  2 and  a corresponding  schematic  is  shown  in  Fig.  3.  The 
driven  burner  tube  with  the  propellant  sample-holder,  the  pressure 
transducer,  and  the  acoustic  driver  are  all  contained  within  the  0 to 
500  psig  pressurization  tank  which  is  equipped  with  high  pressure  hinged 
ports  to  allow  easy  access  for  propellant  sample  changes,  burner  tube 
removal,  and  maintenance. 

Air  supply  for  tank  pressurization  and  the  flow  requirements  of 
the  electropneumatic  driver  is  provided  by  a 3000  psig,  500  cu.ft. 
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Figure  2.  Photograph  of  Pressurized  Driven  Burner  Facility. 


blow  down  facility.  The  tank  pressure  and  driver  airflow  are  main- 
tained by  a pressure  control  valve  in  the  exhaust  line  of  the  system 
as  shown  in  the  schematic  of  the  flow  system  in  Fig.  4. 

The  components  of  the  high  pressure  driven  burner  tube  facility 
and  the  principles  of  its  operation  are  basically  unchanged  from  the 
low  pressure  facility  previously  in  operation.  The  impedance  tube 
contained  in  the  high  pressure  tank  was  fabricated  from  a stainless 
steel  pipe  with  a two  inch  inside  diameter.  Provisions  for  a propellant 
sample  holder  are  included  on  one  end  and  an  acoustic  driver  capable 
of  developing  4,000  watts  of  acoustic  power  is  close-coupled  to  the 
tube  on  the  other  end.  Provisions  for  instrumentation  (pressure  and/or 
temperature  measurements)  have  been  included  along  six  feet  of  the 
tube  wall.  The  instrumentation  locations  are  one-half  inch  apart 
measured  from  the  face  of  the  propellant  sample.  A sketch  of  the 
driven  burner  tube  and  details  of  the  propellant  sample  holder  are 
presented  in  Figs.  5 and  6 respectively.  During  a test,  the  propel- 
lant sample  is  ignited  by  a nichrcme  wire  glued  to  the  sample  surface 
in  an  "S"  shape.  A photograph  of  a propellant  sample  in  a sample  holder 
with  the  nichrome  wire  glued  in  place  is  shown  in  Fig.  7.  To  assure  a 
rapid  and  uniform  ignition  of  the  propellant,  the  surface  of  the  sample 
is  coated  with  a thin  pyrotechnic  mixture  comprised  of  72.4%  potassium 
perchlorate ,14.8%  titanium,  6.9%  Boron  powder,  and  6.0%  polyisobuty lene 
binder  dissolved  in  toluene.  For  ignition  the  nichrome  wire  is  heated 
by  a 40  volt,  15  ampere  power  supply.  Measuring  the  standing  wave  pat- 
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measurement  at  high  temperatures  poses  the  problem  of  protecting  the 
highly  sensitive  transducer  diaphragms  from  the  hot  gases  in  the  tube. 
Various  methods  for  the  protection  of  the  transducers  have  been  in- 
vestigated, The  semi-infinite  tube  technique  described  in  Ref,  2 has 
proven  unsatisfactory  for  high  pressure  operation  because  of  low  fre- 
quency resonance  in  the  infinite  tube  at  high  pressures.  A sketch  of 
an  alternate  pressure  transducer  installation  utilizing  a short  adap- 
tion is  presented  in  Fig,  8,  This  adaption  has  a 3/16-inch  inside 
diameter  and  it  measures  5/8 -inch  from  the  transducer  face  to  the 
inside  of  the  tube  wall.  Further  protection  for  the  transducer  is 
provided  by  a thin  coating  of  clear  silicone  grease  applied  to  the 
transducer  diaphragm.  Extensive  calibration  and  testing  of  this 
short  transducer  adaptor  has  been  conducted  and  the  results  indi- 
cate that  in  the  frequency  range  of  interest  for  this  experiment, 

(i.e.,  300  Hz  to  1100  Hz),  the  adapter  provides  highly  satisfactory 
pressure  amplitude  and  phase  measurement  characteristics.  A detailed 
discussion  with  empirical  data  and  design  criteria  for  short  probes 
to  be  used  for  dynamic  measurements  is  outlined  in  Refs.  7 and  8. 

The  acoustic  driver  used  in  the  high  pressure  driven  burner  tube 
is  a Ling  ETP-94  B electropneumatic  driver  capable  of  developing  4,000 
watts  of  acoustic  power  which  provides  a sound  pressure  level  of  170 
decibels  in  the  burner  tubes.  The  driver  is  close-coupled  to  the  tube 
waLl  at  the  downstream  end  of  the  burner  tube.  The  frequency  and  wave- 
form output  of  the  acoustic  driver  is  controlled  by  a Spectral  Dynamics 
Oscillator,  Model  SD104A-5.  For  this  investigation,  the  waveforms  of  the 
pressure  oscillations  are  sinusoidal  and  the  frequency  of  the  oscillations 
are  constant  during  a given  test. 
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B.  Test  Procedures 


The  data  acquisition  system  and  the  equipment  to  establish  and 
monitor  test  conditions  are  located  in  a laboratory  area  adjacent 
to  the  room  housing  the  high  pressure  driven  burner  tube  experiment. 

A photograph  of  the  control  area  for  the  experiment  is  presented  in 
Fig.  9.  Test  conditions  are  established  by  slowly  pressurizing  the 
high  pressure  tank  housing  and  the  burner  tube  to  the  desired  oper- 
ating pressure.  The  pressurization  tank  pressure  and  driver  airflow 
requirements  are  maintained  by  the  tank  pressurization  control  valve 
in  the  exhaust  line.  With  tank  pressure  stabilized  and  a standing 
waveform  of  a desired  frequency  established  a test  run  is  initiated. 

The  data  acquisition  period  of  a test  run  includes  four  phases;  a 
brief  pre-ignition  test  period  with  the  acoustic  drivers  on  and  test 
conditions  established  in  the  burner  tube,  ignition  of  the  propellant 
sample,  the  propellant  "quasi-steady"  burning  period  and  the  propel- 
lant extinguishment  phase.  The  data  acquisition  period  is  normally 
about  two  and  one-half  seconds.  A more  detailed  discussion  of  the 
data  acquisition  and  data  reduction  procedures  are  presented  in  the 
following  section. 

C.  Minicomputer-Based  Data  Acquisition  System 

During  the  reporting  period  a minicomputer-based  data  acquisition 
system  was  incorporated  in  the  program  to  obtain  improved  pressure  and 
temperature  data.  This  system  processes  the  data  in  three  stages.  First, 
during  a test  the  analog  signals  from  the  transducer  and  thermocouple 
channels  are  sampled,  digitized,  and  stored  at  a controlled  rate.  Next, 
the  stored  readings  are  processed  by  the  computer.  The  transducer  data 
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are  digitally  filtered  and  pressure  amplitude  and  phase  data  are 
obtained  at  the  frequency  of  interest.  Likewise,  thermocouple  signals 
are  converted  to  temperatures.  Finally,  the  pressure  and  temperature 
data  are  printed  out  and  plotted.  These  values  are  then  used  in  the 
data  reduction  scheme  described  in  Section  II  to  obtain  the  admit- 
tance values.  In  the  remainder  of  this  section,  a more  detailed  de- 
scription of  the  data  acquisition  system  is  presented. 

Phase  1:  Data  Sampling.  As  shown  in  Fig.  10  the  data  sampling  phase 
involves  four  major  components  - a Preston  GMAD-11  analog-to-digital 
converter  (ADC),  an  HP  3220-A  frequency  synthesizer  (pacer),  an  HP 
2100S  minicomputer,  and  a disc  storage  unit.  The  ADC  is  capable  of 
sampling  up  to  12  channels  of  input  data  by  use  of  a multiplexer 
and  has  a range  of  + 2.56  volts  with  a resolution  of  20  mv.  To  ob- 
tain maximum  signal  resolution,  the  input  from  the  transducer  and 
thermocouple  channels  are  amplified  so  that  the  maximum  anticipated 
amplitude  for  a test  is  as  close  to  2.56  volts  as  possible  without 
exceeding  this  upper  limit. 

Upon  initiation  of  a test  the  following  sequence  of  events  occurs 
as  denoted  in  Fig.  10.  (1)  The  computer  initiates  a signal  to  start 
the  analog-to-digital  conversion  process.  (2)  The  pacer,  capable  of 
controlling  the  sampling  rate  up  to  12  million  samples /second,  also 
initiates  a signal  at  the  desired  sampling  rate.  In  the  present  ex- 
periment up  to  12,000  samples  per  second  are  taken.  When  the  pacer 
and  computer  signals  are  coincident  the  sample  and  hold  channels  of 
the  ADC  go  into  the  "hold"  condition  so  that  the  input  signal  voltages 
at  that  instant  of  time  are  retained.  (3)  The  starting  signal  from  the 


switch  proceeds  to  the  multiplexer  which  indexes  through  the  sample- 
and-hold  channels.  (4)  If  the  data  from  all  the  specified  channels 
have  been  digitized  (determined  by  the  divide  logic  circuitry),  then 
more  data  from  the  input  channels  are  sampled  and  held.  (5)  The  data 
from  the  sample-and-hold  device  are  sequentially  transfered  two 
channels  at  a time  by  the  multiplexer.  (6)  The  voltages  from  the 
first  of  the  two  channels  sampled  by  the  multiplexer  enters  an  8-bit 
voltage  to  binary  converter,  is  digitized,  and  is  stored  in  a 16-bit 
output  register.  (7)  The  second  channel  is  also  digitized  and  stored 
in  the  register.  (8)  When  this  process  is  completed,  the  computer 
stores  the  two  8-bit  bytes  as  a one  16-bit  word.  By  packing  the  data 
in  this  manner  the  data  throughput  rate  is  essentially  doubled  and 
increases  the  upper  frequency  limit  at  which  reliable  data  can  be 
acquired.  With  the  present  system,  reliable  pressure  data  at  fre- 
quencies of  up  to  6000  Hz  can  be  obtained.  Potentially,  the  upper 
frequency  limit  can  be  extended  to  one  mega  Hz.  (9)  The  packed  data 
word  is  stored  in  a buffer  in  the  computer,  and  the  sequence  of  events 
is  repeated  starting  with  (1). 

Before  testing,  the  sampling  rate,  the  size  of  the  data  buffer, 
and  the  sampling  interval  are  determined.  The  sampling  rate  depends 
upon  the  frequency  of  the  driven  oscillations.  The  samples  are  taken 
rapidly  enough  to  obtain  a sufficient  number  of  data  points  per  cycle 
period  for  good  signal  definition.  The  data  buffer  and  sampling  inter- 
val are  determined  by  the  test  conditions.  In  this  study,  the  ignition 
and  burnout  transients  produce  large  variations  in  pressure  amplitude 
and  phase.  The  sampling  interval  must  be  small  enough  to  detect  these 
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variations  so  that  they  can  be  distinguished  from  the  nearly  constant 
amplitudes  which  occur  during  the  period  of  quasi-steady  burning. 

On  the  other  hand  the  sanpling  interval  must  be  long  enough  to  ensure 
that  the  digital  filtering  interval  is  sufficient.  This  last  require- 
ment can  be  relaxed  for  the  tests  conducted  in  the  present  investi- 
gation since  the  signals  thus  far  encountered  have  been  relatively 
free  of  background  noise.  Once  the  sampling  interval  is  established 
the  size  of  the  buffer  in  the  computer  into  which  the  data  is  stored 
is  determined. 

During  a test  data  are  taken  at  up  to  80  sampling  intervals  of 
from  10  to  50  milliseconds,  and  the  time  between  intervals  is  from 
15  to  30  milliseconds.  The  duration  of  the  run  is  from  two  to  six 
seconds  for  the  3/8"  -thick  propellant  samples  used  and  from  20  to 
40  data  points  are  obtained  during  steady  burning.  After  the  digi- 
tized data  points  taken  during  one  sampling  interval  have  been  stored, 
they  are  then  written  onto  a disc  so  that  room  can  be  made  in  the 
buffer  to  receive  the  data  from  the  next  sampling  interval.  These 
steps  are  depicted  in  Fig.  11. 

Phase  2:  Data  Processing.  During  the  data  processing  phase,  the  disc, 
computer,  and  a Tektronix  4012  graphics  display  terminal  are  used. 

The  steps  involved  are  shown  in  Fig.  12.  (1)  The  digitized  data  taken 
over  one  time  interval  is  read  from  the  disc  into  the  computer.  (2) 

The  computer  software  unpacks  the  data  words  so  that  the  data  from  each 
channel  can  be  processed.  From  the  digitized  signals  of  the  transducer 
channels,  the  average  amplitude  and  phase  over  the  time  interval  is  ob- 
tained. If  the  time  interval  is  sufficiently  short,  the  time  variations 
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Figure  11.  Signal  Processing  During  Data  Sampling  Phas 
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in  amplitude  and  phase  will  not  be  significant  and  the  computed 
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averages  will  be  close  to  the  true  value.  The  average  temperature 
over  the  time  interval  is  also  computed  from  the  digitized  thermo- 
couple data  by  converting  the  measured  voltages  to  temperature  using 
a Forsyth  polynomial  expansion  obtained  from  the  thermocouple  cali- 
bration charts.  Steps  one  and  two  are  repeated  for  each  successive 
time  interval.  Step  (3)  consists  of  taking  the  resulting  amplitude 
and  phase  data  and  plotting  them  against  nondimensional  time  on  the 
terminal.  This  provides  a check  to  ensure  that  proper  test  conditions 
were  established  during  the  run.  The  data  processing  phase  generally 
requires  less  than  five  minutes  of  execution  time  so  that  if  desired 
test  conditions  were  not  achieved,  minimal  time  is  lost  between  the 
preliminary  data  check  and  retesting.  In  addition  to  being  displayed 
on  the  terminal,  the  amplitude,  phase,  and  temperature  data  are  also 
stored  on  one  disc  track  preassigned  to  the  test  which  is  protected 
by  software.  These  data  can  then  be  displayed  at  a later  date  so  that 
minimal  time  between  runs  is  lost. 

Phase  3:  Data  Presentation.  The  data  presentation  phase  simply  consists 
of  reading  the  data  stored  during  a run  off  the  disc  track  assigned  to 
that  test  and  printing  or  plotting  the  results.  The  first  step  is  to 
plot  the  amplitude  and  phase  against  nondimensional  time.  From  these 
plots  the  time  intervals  at  which  steady  state  burning  was  achieved 
can  be  ascertained.  At  these  selected  points  the  amplitude,  phase,  and 
temperature  data  can  - be  plotted  or  printed  out  versus  distance. 
These  data  are  then  used  in  the  analyses  described  in  Section  II  to 
compute  the  admittance.  This  final  admittance  computation  is  carried 
out  on  a CYBER-70  computer  instead  of  the  minicomputer  because  of  the 
time  and  storage  requirements  involved. 
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IV.  EXPERIMENTAL  RESULTS  AND  SUMMARY 


The  experimental  efforts  during  this  reporting  period  were  directed 


toward  improving  the  experimental  techniques  and  determining  the 


admittances  of  a number  of  solid  propellants  at  high  pressure  test 


conditions.  To  date,  twenty  nine  tests  have  been  conducted  using 


three  different  propellants:  (1)  a T-13  propellant,  (2)  an  undesig 


nated  white  propellant  supplied  by  Thiokol  (referred  to  as  propellant 


B in  this  report),  and  (3)  an  A-15  propellant.  A summary  of  the  tests 


conducted  thus  far  is  presented  in  Tables  II  through  IV 


TABLE  II.  Tests  Conducted  with  Propellant  B 


Propellant 

Type 


Driving 
Frequency 
in  Hz 

Chamber 
Mean  Pressure 
in  psig 

625 

300 

278 

300 

502 

300 

742 

300 

967 

300 

742 

300 

967 

300 

742 

300 

985 

300 

1 

Closed 

I 

Open 

II 

Open 

II 

Open 

II 

Open 

II 

Closed 

I 

Closed 

I 

Muffler 

End 

IV 

Muffler 

End 

IV  1 

TABLE  III.  Tests  Conducted  with  T-13 


Expt.  No. 
E-16-633 

Propellant 

Type 

Driving 
Frequency 
in  Hz 

Chamber 
Mean  Pressure 
in  psig 

Exhaust  End 
Configuration 

551 

T-13 

625 

300 

Closed 

I 

560 

II 

278 

300 

Open 

II 

562 

II 

502 

300 

Open 

II 

564 

ft 

742 

300 

Open 

II 

566 

" 



967 

300 

Open 

II 

TABLE  IV.  Tests  Conducted  with  A-15 


E-16-633 


Driving 
Frequency 
in  Hz 

Chamber 
Mean  Pressure 
in  psig 

625 

300 

967 

300 

967 

300 

1128 

300 

967 

300 

300 

1128 

300 

980 

300 

912 

300 

181 

300 

181 

300 

967 

0 

181 

0 

181 

0 

967 

0 

Exhaust  End 
Configuration 


Closed  I 

Closed  I 

Closed  I 

Closed  I 

Open  II 

Open  II 

Open  II 

Extended  Open  III 

Extended  Open  III 

Closed  I 

Muffler  End  IV 

Open  II 

Open  II 

Closed  II 

Closed  tt 


A schematic  of  the  exhaust  end  configurations  used  during  the 
tests  are  presented  in  Fig.  13.  Different  configurations  were  em- 
ployed to  ensure  that  the  conditions  at  the  exhaust  end  did  not 
affect  the  relative  pressure  amplitudes  and  phases  measured  along 
the  standing  wave,  from  which  the  admittance  of  the  propellant  is 
obtained.  The  data  obtained  to  date  indicate  that,  to  within  ex- 
perimental error,  no  significant  change  in  the  admittance  values 
occurs  with  changing  exhaust  end  conditions. 

In  Figs.  14  through  17,  typical  data  is  presented  for  a test 
conducted  at  742  Hz  with  propellant  B with  an  open-ended  exhaust 

condition.  Figure  14  is  a plot  of  the  pressure  amplitude  in  decibels 

-4 

(re.  2x10  p bar)  versus  record  number  for  four  pressure  channels 
located  at  different  axial  stations  along  the  tube.  Each  record 
represents  a short  time  interval  during  which  data  was  recorded 
by  the  analog-to-digital  converter.  Consecutive  records  correspond 
to  consecutive  time  instances;  hence  the  abscissa  can  be  considered 
as  representing  time. 

The  data  in  Fig.  14  can  be  divided  into  four  phases:  (1)  the 
pre-ignition  period  during  which  the  Ling  acoustic  drivers  are  the 
only  source  of  wave  excitation  (records  1 to  14);  (2)  the  propellant 
ignition  phase  (records  15  to  18);  (3)  the  quasi-steady  burning  phase 
during  which  reliable  pressure  data  can  be  taken  to  obtain  the  burning 
solid  propellant  admittance  values  (records  19  to  40);  and  the  burnout 
phase  (records  41  to  80).  During  the  ignition  and  burnout  phases,  the 
temperature  in  the  tube  is  changing  rapidly  with  time  which  produces 
large  variations  in  the  wavelength  of  the  standing  wave.  Therefore, 
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reliable  data  cannot  be  obtained  during  these  periods.  On  the  other 
hand,  no  significant  rapid  temperature  variations  are  expected  during 
the  quasi-steady  phase  of  the  test.  Hence,  data  measured  during  this 
period  should  produce  reliable  admittance  values  for  the  burning 
solid  propellant. 

The  behavior  of  the  phase  with  respect  to  time  is  shown  in  Fig.  15. 
All  phases  shown  in  this  plot  are  referenced  to  the  phase  measured 
three  inches  from  the  propellant  surface.  If  the  phase  differences 
are  positive,  then  the  pressure  oscillations  seven,  nine  and  ten 
inches  from  the  sample  lead  the  referenced  pressure.  This  type  of 
condition  indicates  that  the  incident  wave  initiated  by  the  acoustic 
driver  is  greater  in  amplitude  than  the  wave  reflected  from  the  burn- 
ing propellant  surface;  therefore,  the  combustion  process  is  damping. 

On  the  other  hand,  if  the  pressure  oscillations  further  from  the 
burning  propellant  lag  the  oscillations  closest  to  the  sample,  then 
the  phase  difference  is  negative  and  the  wave  reflected  from  the  pro- 
pellant is  greater  in  amplitude  than  the  wave  incident  from  the  acoustic 
driver.  This  condition  indicates  that  the  combustion  process  is  driving. 

From  the  phase  data  in  Fig.  15  it  appears  that  the  combustion  pro- 


cess is  randomly  damping  and  driving  the  acoustic  oscillations  during 
the  run.  This  result  is  unexpected  and  behavior  of  this  type  has  not 
been  reported  in  other  studies  of  burning  solid  propellants  surveyed 
under  this  contract.  The  pressure  amplitude  and  phase  variation  with 


distance,  taken  at  record  20,  are  presented  in  Figs.  16  and  17.  The 


■ «-• 


data,  an  admittance  value  of  0.1  was  computed  which  means  that  the 
combustion  process  increased  the  incident  wave  amplitude  by  approxi- 
mately twenty  per  cent  which  indicates  a strong  driving  condition 
at  the  propellant  surface.  However,  large  damping  conditions  are  ob- 
served at  other  times  during  the  burn  (e.g.  see  Record  35).  An  investi- 
gation of  the  phenomenon  is  currently  underway. 

In  Figs.  18  through  21  the  behavior  of  the  amplitude  and  phase 
of  the  pressure  oscillations  with  respect  to  time  and  distance  are 
presented  for  the  T-13  propellant  at  a frequency  of  742  Hz  with  an 
open-end  exhaust  condition.  The  two  "spikes"  in  the  data  in  Figs.  18 
and  21  at  records  38  and  58  are  erroneous  data  caused  by  the  ADC 
digitization  process  when  a significant  bit  is  either  erroneously 
included  or  excluded.  As  seen  from  Fig,  18,  the  pressure  amplitude 
rises  and  falls  periodically  with  time  during  the  burn.  Since  this 
behavior  occurs  only  after  ignition,  it  is  apparently  associated 
with  the  combustion  process.  From  the  phase  data  presented  in  Fig. 

19,  no  large  driving  or  damping  characteristics  are  exhibited  by 
the  combustion  of  the  T-13  propellant.  Instead,  it  behaves  essen- 
tially like  a hard  termination  and  does  not  increase  or  decrease  the 


amplitude  of  the  incident  wave  significantly.  The  axial  variation  of 
the  pressure  amplitude  and  phase  at  record  45  is  shown  in  Figs.  20 
and  21.  From  the  phase  plot.  Fig.  21,  a weak  damping  by  the  combus- 


tion process  is  exhibited.  However,  at  other  times  (i.e.  record  50) 
a weak  driving  condition  is  shown. 

In  almost  all  tests  conducted  thus  far,  at  300  psig,  the  tested 


propellants  randomly  exhibited  driving  at  some  times  and  damping  at 
others.  The  propellants  which  are  the  most  unstable  (i.e.,  the  A-15 
and  propellant  B)  seem  to  exhibit  the  largest  variation  from  driving 
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to  damping  whereas  the  most  stable  propellant  (i.e.  the  T-13)  seems 
to  show  the  least  variation.  Current  research  efforts  are  concerned 
with  the  development  of  an  appropriate  explanation  for  the  measured 
experimental  data. 
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Data  taken  at  300  psig  indicate  that  the  combustion  process  of  the  o.ic 
lant  periodically  drives  and  damps  acoustic  oscillations  under  most  of  the 
conditions  encountered. 


